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ABSTRACT 

Using direct numerical simulations, we demonstrate that "flow IV" of Roberts (1972) ex- 
hibits dynamo action dominated by horizontally averaged large-scale magnetic field. With 
the test-field method we show that the magnetic eddy diffusivity is negative and overcomes 
the molecular diffusivity to explain this large-scale dynamo for magnetic Reynolds numbers 
above sa 8. As expected for a dynamo of this type, the two horizontal field components grow 
independently of each other and have arbitrary amplitude ratios and phase differences. Small 
length scales of the mean magnetic field are shown to be stabilized by the eddy diffusivity be- 
coming positive at wavenumbers larger than twice the characteristic wavenumber of the flow. 
For magnetic Reynolds numbers below « 0.5 the eddy diffusivity is confirmed to be positive, 
as expected for all incompressible flows. Earlier claims of a dynamo driven by a modified 
Taylor-Green flow through negative eddy diffusivity could not be confirmed. 
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1 INTRODUCTION 

The equations of magnetohydrodynamics (MHD) permits 
growth of magnetic energy at the expense of kinetic en- 
ergy. This phenomenon is called the dynamo effect (see e.g. 
[Brandenburg & Subramaniaiil l2005l for a recent review). If the 
dynamo effect gives rise to a magnetic field whose characteristic 
length scale is greater than that of the fluid, we call it a large-scale 
dynamo. Most astrophysical dynamos, including the solar dynamo 
and the galactic dynamo, are of this type. 

To theoretically describe the large-scale dynamo one must av- 
erage the equations of MHD over the small scales to write an ef- 
fective equation for the large-scale magnetic field. This effective 
equation can be wri tten down by using either mean-field theory 
jSteenbeck, Krause'& Radler 1966) or multiple-scale expansions 
(see e.g. Zheli govskvl I20 1 ll , for a recent review). These equations 
contain turbulent transport coefficients: the a effect and turbulent 
(or eddy) diffusivity. In general, both are tensors whose complexity 
depends on the symmetries of the problem. Within the formalism of 
mean-field theory, it is generally a non-trivial task to calculate the 
turbulent transport coefficients even if we ignore the back-reaction 
of the magnetic field on the flow, i.e., for kinematic dynamos. 
For several kinematic problems, the turbulent transp ort coefficients 
have been calculated using the test-field method (SuretaL I l2008l : 
[Brandenburg et al1l2008a ,b). Typically, it is found that the a ef- 
fect gives rise to the growth of a large-scale magnetic field while 
the eddy diffusivity contributes to decay by effectively enhancing 
the molecular magnetic diffusivity. However, multiscale methods 
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have shown that for certain flows the a effect can b e zero, but the 
eddy diffusivity may turn out to be negative (see e.g.,'Lanotte et alj 
1999; Zheligovsky et al. 2001). If, for such flows, the magnitude of 
the eddy diffusivity is greater than the molecular one, a large-scale 
dynamo can operate. 

In this paper we first demonstrate, by direct num erical simu- 
lations (DNS) that for a particular flow ('Ro berts|[T972l) , which will 
in the following be refeiTed to as Roberts-IV flow, it is possible to 
drive a kinematic large-scale dynamo although the a effect (and the 
planar-averaged kinetic helicity) is indeed zero. Next, by using the 
test-field method, we show that such a dynamo can be accurately 
described by zero a effect but negative eddy diffusivity which dom- 
inates over molecular one. 



2 THE ROBERTS-IV FLOW 

In a remarkable paper, iRobertsI (Il972h studied several two- 
dimensional periodic flow patterns that all produce dynamo action, 
which means that an infinitesimally weak magnetic field gets am- 
plified expo nentially in time. In connection with understanding the 
geodynamo, lTilgneil|2004.) studied in some detail the "flow IV" of 
lRobertsl ( fT973) . We follow here Tilsner's definition of the flow: 



U 



J2/ f sin kox cos koy 

"^0 I —yj2/f cos fcoa; sin fcoy 
sin kox 



(1) 



where uo characterizes the amplitude of the flow. Its vorticity, 
W — 'V X U,is given by 
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Figure 1. Three components of the magnetic field on the periphery of the 
computational domain for Rm = 20, / = 1, and domain size Lx = Ly = 
Lz = 27r/fco. 
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= V X {U X B -T^J), 



(4) 



yT" ko cos kox 
2/f ko sin kox sin koy ) 

Here, the parameter / determines tiie relative importance of eddy 
diffusion along the z direction (large /) compared to that in the xy 
plane (small /). The kinetic helicity density, W ■ U, is given by 

W ■ U — V2uQko (1 + sin'^fcox) sinfcoy (3) 

and is independent of /. lTilgneilll2004l) showed that in spite of the 
horizontally averaged kinetic helicity density, W ■ U, being zero; 
the Roberts IV flow gives rise to dynamo action. In other words, it 
leads to growing solutions of the induction equation, 

dB 

dt 

where rj is the microphysical (molecular) magnetic diffusivity, B 
is the magnetic field, J = V x B is the current density, and we 
have chosen our units such that the vacuum permeability is unity. 
Here and in the followin g, overbars den ote xy averaging. 

Note, however, that lTilgneJ (|2004|) described the dynamo to be 
a small-scale one, i.e., the characteristic length scales of the mag- 
netic field is of the same order as 1/fco . In the following, we obtain 
solutions to Equation © via DNS using the PENCIL COD^. We 
do not evolve the flow, hence we study kinematic dynamo solutions. 
As an example of the resulting magnetic field, we show in Figure[T] 
the three components of the magnetic field at the periphery of the 
computational domain. It is remarkable that the resulting magnetic 
field has a large-scale component that survives horizontal averag- 
ing, i.e., B = B{z,t) is non-vanishing; see Figure |2] where we 
show examples of the resulting mean field obtained by averaging 
the solution of the DNS. In other words, we have here an example 
of a mean-field dynamo, where i? is a solution of the horizontally 
averaged induction equation. 



dB 

'dt 



= V X (U X B + £ --nJ) 



(5) 



where £ — u x b is the mean electromotive force resulting from 
correlations of fluctuating velocity and magnetic fields, u — U — 
U and b = B B, respectively. Empirically, we find that the 
horizontally averaged solutions of Equation (|4) are of the form 

(Bqx cos{kz + (f>x)\ 
Boy cos{kz + (j)y) e^*, (6) 
/ 

where Bqx, Boy, (px, and (j)y are arbitrary constants, i.e., the x and 
y components of the magnetic field evolve independently of each 
other and they have arbitrary phase shifts, depending just on the 
properties of the initial conditions; see Figure |2] for an example. 
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Figure 2. Examples of runs with different initial conditions showing the 
X and y components of the mean field (normalized by the rms value of the 
actual field) versus z, obtained from the DNS for / = 1, Lx = Ly = 
Lz = 27r/fco, and rj = 0.05 «o /'so, con'esponding to Rm = 20. 



Solutions of Equation ^ can be obtained by mean-field simula- 
tions (MFS) which requires a closed expression for £ in terms of 
U and B. This will be discussed in th e following. 

As pointed out by iTilgneJ ( 12004) ■ the Roberts-IV flow has 
no a effect. This is also confirmed by our simulations. As this 
is a laminar flow, driving a dynamo via fluctuations of a effect 
(Mitra & Brandenburg 2012) is also not possible. This suggests 
that the observed mean fie ld might be produced by a negative e ddy 
diffusivity effect I lLanotte et aljri999l : IZheIigovskv et al.ll200lh . To 
investigate such a possibili t y, we n ow ap ply the test -field method 
(TFM) of Schrinner et al. (2005, 200T); see also 'BrandenburgI 
(2005); Brandenburg et al. (2008a, b) for its implementation in pe- 
riodic geometry; to calculate the turbulent transport coefficients of 
the Roberts IV flow. 

As is long recognized ( lRadleilll976h . the connection between 
£ and i? is a nonlocal one that is described by a convolution of the 
form 



£■1 = Q-i] ° Bj — fjij O Jj, 

where "o" denotes a convolution in space and time, i.e.. 



fj.j O Jj = 



,t — t') Jj {z' , t') dz' dt' , 



(7) 



(8) 



and like wise for , o , , but this term is vanishing for the Roberts- 
IV flow j Tilgneill2004l) . The hats on Aij and fjij indicate that the 
corresponding quantities are integral kernels. 

We emphasize that in Equation (|7j we have made use of the 
fact that the only non-vanishing derivatives of a horizontally av- 
eraged mean field are dBx/dz and dBy/dz, which and be ex- 
pressed in terms of components of J, so the corresponding eddy 
di ffusion tensor is just of rank 2, not, as in the general case of rank 
3 dKrause & Radlej|l9 80'). 

In the TFM, the kernel formulation above is most naturally 
considered in Fourier space with 



^ |http : //pencil- code . googlecode . com/| 



£{k,uj) — aij{k,u})Bj{k,uj) — fjij{k,u})Jj{k,u}), 



(9) 
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where tildes denote appropriately normalized Fourier tra nsforma- 
tion of the corresponding mean-field q uantities (Branden burg et alj 
l2008bl : lHubbard & Brandenburg|l2009l) . 

In the following we consider a three-dimensional domain of 
size Lx X Ly X Lz. For most cases we choose cubic domains, 
i.e., Lx ~ Ly — Lz — 2n/ko. We are interested in the spe- 
cial case of harmonic solutions of Equation (|5} with fixed vertical 
wavenumber fc = fco of a magnetic field that is growing or de- 
caying exponentially proportional to e^*. So we are interested in 
the case uj — iX. As an approximation, we begin by considering 
the case u) = 0, i.e., we ignore the so-called memory effect; see 
iHubbard & Brandenburd ( 12009.) for illustrating the departure in the 
case of the standard Roberts flow with helicity (also known as the 
Roberts-I flow). 



3 RESULTS 

3.1 Sign change of eddy diffusivity 

As we have mentioned already, all components of atj (k, ui) vanish 
for the Roberts-IV flow. Moreover, T)ij{k,u)) is isotropic, i.e., we 
can write f]ij = ??t<5ij. The resulting values of r;t(fco, 0) are shown 
in Figure|3]as a function of magnetic Reynolds number, 

-Rm = uo/rjko. (10) 

For comparison with earlier work, we note that this definition of 
Rm is close to a definition in terms of the rms velocity of the flow 
(for / = 1 we have Umis ~ 1.225 iio) and the wavenumber of 
the energy-carrying eddies kf, i.e., Mrms/7?fcf. If we approximate 
k{ ~ Wrms/wrms, whcfe Wrms is the rms value of the fluctuating 
part of the vorticity, then we have k{ « 1.29 ko- Therefore, we have 
UTvaa/rjki « 0.95 -Rm, which is close to _Rni. 

As is commo n to man y turbulent transport coefficients 
JSur et al.. 2008: Br andenburg et al. 2008b), rjt grows linearly with 
Rm for Rra ^ 0.5; see Figure[3] More importantly, r]t is positive, 
which is to be expected based on a calculation for incompressible 
flows using the second order correlation approximation (SOCA), 
which is valid for Rm «C 1. To show this, one uses the fact that 
the Fourier transform of the velocity correlation tensor is positive 
semidefinite. We note in passing that this is not true for potential 
flows, for which t he eddy diffusion ten sor is negative semidefinite 
for Rm <C 1; see lRadler et al.l ( 1201 ih for a recent demonstration 
using the TFM. 

Returning now to the Roberts-IV flow, which is indeed incom- 
pressible, we show in Figure |3]that i]t changes sign from positive 
to negative at Rm ~ 4. This is clearly a result that cannot be re- 
covered by SOCA. Moreover, for Rm ~ 8, the total magnetic dif- 
fusivity, rj + rit{k, 0) becomes negative, i.e., dynamo action by the 
negative magnetic diffusivity effect is possible. The critical value 
of -Rm agrees with that found above through DNS. The growth rate 
of the dynamo is given in implicit form as a solution of the equation 

X{k) = -[ri + rit{k,iX)]k^ (11) 

for k = fcf). However, it is common to approximate r;t(fco,(^) by 
?7t(fco, 0), and we refer to the corresponding solution as 

Xik) ^ Adisp (fc) = -~[v + Vtik, 0)]k\ (12) 

which is shown in the third panel of Figure|3]and compared with the 
growth rate Adns obtained by solving Equation Q through DNS. 
The agreement between Adisp and Adn s is moderate and remi- 
niscent of what has been found earlier dHubbard & Brandenburd 
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Figure 3. Turbulent magnetic diffusivity versu.s Rm for / = 1 and fc = fco . 

By using a test field that grows exponentially at a rate 
that is equal to the expected growth rate, A — Adns, we find 
for 1] = 0.020 uo/fco (corresponding to Rm ~ 50), the value 
77t(fco,iA) = —0.081 ito/fco with A — 0.061 iiofco, instead of the 
value —0.070 Mo/feo for A = 0. Thus, perfect agreement between 
DNS and TFM is obtained once the memory effect included. 

3.2 Dependence on / 

Let us now discuss the dependence on the parameter /, which char- 
acterizes the relative importance of vertical to horizontal eddy dif- 
fusion. We consider here the case of Rm = 20 and k — ko- The 
results are shown in Figured The negative eddy diffusivity dynamo 
is found to be operating in the range 0.6 < / < 1.23, i.e., when 
the vertical eddy diffusivity is not much larger than the horizontal. 

As indicated by the dotted line in Figured both for small and 
for large values of /, there is an approximate power law depen- 
dence with |j7t| ~ f^'''- However, in the range 1.3 < / < 3 the 
TFM diverges and is unable to deliver useful results. Diverging re- 
sults of the TFM are common and related to unstable eigenvalues 
of the associated homogeneous system of equations solved in the 
TFM. Usually, this problem can be avoided by focussing on the 
inhomogeneo us term, which dominates the resu lt at intermediate 
time intervals dRheinhardt & Brandenburj|2010l) . 
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Figure 4. Dependence of eddy diffusivity on / for Rm = 20 and k = kg. 
Negative values of r]t are indicated with open symbols, and the horizontal 
dash-dotted line indicates the region above which there is dynamo action 
when r; + »7t is negative. The dotted line has a slope of 1.7 and is shown for 
orientation. 
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Figure 5. Wavenumber dependence of eddy diffusivity for Rm = 20 and 
/ = !■ 



3.3 Scale dependence 

Owing to the fc^ factor in Equations JIU and ( I12t . one might ex- 
pect dynamos driven by negative eddy diffusivity to grow faster at 
smaller scales, unless r]t becomes positive at larger values of k. In 
Figure |5] we show the k dependence of rit(k,0). It turns out that 
r;t(fc,0) is approximately constant for k < k* ^ 1.125 fco, and 
positive with an approximate dependence 



r?t(fc,0) 



0.31 uo/ko 



(for k — fcoo), 



(13) 



where koo ~ 2.25 fco. In the range k, < k < k^o we have two data 
points in Figure |5] that clearly deviate from the description above. 
In addition, there are several other values of k in this range where 
the TFM again diverges and is unable to deliver useful results. 

To illuminate the problem of intermediate k values further, we 
now use DNS to compute the growth rate as a function of the do- 
main size Lz, keeping Lx = Lz = 27r/fco. The result is shown in 
Figure|6] It turns out that A {— Adns) has a maximum at Lz ko ~ 3. 
Thu s, large scale separ ation, as assumed in some analytic stud- 
ies jLanotteet all 19991) , is neither needed nor necessarily helpful 



Figure 6. Dependence of the growth rate A and the oscillation frequency ui 
on the domain size Lz for i?m = 20 and / = 1. In all cases with A > 0, 
we have uj = 0. 



for the operation of this negative eddy diffusivity dynamo. Further- 
more, for Lzko < 2.55, no dynamo is possible and the field decays 
in an oscillatory fashion. The oscillation frequency uj grows sharply 
as Lz decreases further; see the dotted line in Figure |6] 



4 MODIFIED TAYLOR-GREEN FLOWS 

While the possibility of a dynamo driven by negative eddy diffu- 
sivity has not been previously quantified f or the R oberts-IV flow, 
it was discussed in some detail by .Lanotte et alj 17999,) for the 
Taylor-Green (TG) flow, 

(sinfcoxcosfcoj/cosfcoz \ 
— cos fcoa; sin fco J/ cos fco 2 1 , (14) 
/ 

and the modified TG flow, Utg + AUa + BUb, where A and B 
denote the amplitudes of additional contributions proportional to 

(sin 2fco3; cos 2fco z \ 
sin 2fco2/ cos 2fco2: 1 
— (cos 2fcoa; + cos 2fcot/) sin 2fco2 / 



and 



(sin fcoa; cos 3koy + A sin 3fcoa; cos fco y) cos koz 
-{cosSkoxsinkoy + ^ cosfcoxsin3fcoy) cosfcoz 
(cos fcoa; cos 3koy — cos 3fco2; cos koy) sin koz 



respectively. 

We perform kinematic DNS with this flow and we indeed find 
dynamo action. We carry out calculations for ri = 0.02uo/fco, 
which corresponds to the case where the dynamo is mildly super- 
critical. We consider the following three cases: (a) A cube of size 
Lx = Ly = Lz = 27r/fco, using 128"^ meshpoints; (b) a cuboid 
with Lz = 4Lx and Ly — Lx = 27r/fco, using 128^ x 512 mesh- 
points; and (c) a cuboid with Lx ~ Ly — ALz and Lz = 2n/ko. 
In all these three cases the large-scale field obtained by averag- 
ing over the xy plane decays as a function of time, i.e., no large- 
scale dynamo is obtained. However, there still remains, in principle, 
the possibility of a large-scale field developing that is zero under 
Reynolds averaging but is non-zero under other averaging proce- 
dure, e.g., Fourier filtering. But even this possibility is ruled out 
because we observe no growth of a large-scale field at fc < fco; 
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Figure 7. Power spectra of DNS for the TG flow with A = 1 and B = 
in a domain of size = Ly = = Sir/ko using 256"' mesh points and 
T] = 0.02uo/ko at different times during the exponential growth phase of 
the dynamo. The earliest time is shown as a dash-dotted line. 

see Figure|7] This is found by calculating the spectrum of the mag- 
netic field from simulations with 256"^ meshpoints and a domain 
size of Lx = Ly = Lz = 87r/fco. In other words, a dynamo is 
observed but it is not a large-scale dynamo. This is corroborated by 
the TFM which produces positive eddy diffusivity in all the above- 
mentioned cases. In case (a) rjt ~ +0.135 uo/^o for A — 1 and 
B = 0, and rjt — +0.146 ito/fco for A = 1 and B = 1; in case 
(b) »7t — +0.158 uo/fco. However, in case (c) the TFM becomes 
unstable. These result are in complete contrast with the Robeits-IV 
flow, for which the corresponding values are of the order of unity; 
see Figure |2] 



5 CONCLUSIONS 

In the present work we have revisited the Roberts-IV flow using 
the TFM to compute the full set of turbu l ent tra nsport coefficients. 
We confirm an earlier result of iTilgned (|2004 that a dynamo is 
possible and that all components of the a tensor are vanishing. In 
addition, we find that there is a finite horizontally averaged mean 
magnetic field, which should be explicable in terms of mean-field 
dynamo theory. The TFM reveals that the eddy diffusivity tensor 
is diagonal. Moreover, in the regime where the dynamo is excited, 
the eddy diffusivity is sufficiently strongly negative such that the 
total (molecular plus eddy) diffusivity is negative. This is an un- 
usual situation in that the horizontal components of the mean field 
are completely decoupled and grow independently with arbitrary 
relative amplitudes and phase shifts. 

Many laminar flows are only slow dynamos, and the Roberts- 
IV flow is no exception. These dynamos are therefore not expected 
to be astrophysically relevant. However, the method used to study 
such dynamos (TFM combined with DNS) is now playing an im- 
portant role in the study of astrophysical dynamos for turbulent 
flows. The present work highlights the accuracy of this method in 
that it enables us to pinpoint the detailed nature of a dynamo ex- 
hibiting a finite averaged magnetic field. 

In the present case of laminar flow patterns, nonlocality is 
crucially important. In other words, turbulent transport is de- 
scribed by a convolution with the mean fields in space and time 
rather than just a multiplication. The TFM is particularly well 
suited to deal with such cases. For generic turbulent flows, as 



shown in earlier works by 'Hu bbard & Brandenburg! (|2009|) and 
Rheinhardt & Brandenburg (201^, we expect these transport ker- 
nels to have a relatively simple form and th at complicated kernels, 
such as found here and in the earlier work dRadler & Brandenburg! 
!2009l) are atypical. Note however, that even though most astro- 
physical flows are turbulent and are expected to become statisti- 
cally homogeneous and isotropic at small scales; in practice large 
scale anisotropy and inhomogeneity plays an important role. In 
many of those cases nonlocality cannot be neglected and many 
Fourier modes need to be taken into account, as demonstrated by 
Chatteriee et al.! JioT l|) for flows driven by the magnetic buoyancy 
instability. 
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